function [] = plot_WAcoast(type, varargin)
%-----------------------------
%
% plot_WAcoast(type)
%
% this adds a coastline to the current plot of type:
% 'detailed': loads coastline_detailed.mat which is good for Puget Sound
% 'regional': loads coastline_regional.mat which is good for WA coast and S of Georgia
%
% varargin are line specs to pass to plot 
%
% ex) 'linewidth',2,'color','k' (default is black)
%-----------------------------


%% Set global paths and directories
addpath('../../alpha/');
Tdir = toolstart;
addpath(Tdir.data);


%% Default plot type
if nargin < 1 %default to regional
    type = 'all_noBorder';
end


%% Load the correct coastline and bathymetry files to use
dirname = [Tdir.data 'geo_data/'];
switch type
    case 'detailed'
        filename = [dirname, 'coast/pnw_coast_detailed.mat'];
        load(filename);
    case 'regional'
        filename = [dirname, 'coast/pnw_coast_regional.mat'];
        load(filename);
    case 'all'
        filename = [dirname, 'coast/pnw_coast_combined.mat'];
        load(filename);
        lon_coast = lon;lat_coast = lat;
    case 'all_noBorder'
        filename = [dirname, 'coast/pnw_coast_combined_noborder.mat'];
        load(filename);
        lon_coast = lon;lat_coast = lat;
    case 'all_wbathy'
        filename = [dirname, 'coast/pnw_coast_combined.mat'];
        load(filename);
        lon_coast = lon;lat_coast = lat;
        filename = [dirname, 'hybrid/pnw_combined_full.mat'];
        load(filename);
        contour(lon_topo,lat_topo,z_topo,[-4000 -3000 -2000 -1000 -500 -200 -100 -50],'color',0.6.*[1 1 1])
    case 'regional_wbathy'
        filename = [dirname, 'coast/pnw_coast_regional.mat'];
        load(filename);
        filename = [dirname, 'hybrid/pnw_combined_full.mat'];
        load(filename);
        contour(lon_topo,lat_topo,z_topo,[-4000 -3000 -2000 -1000 -500 -200 -100 -50],'color',0.6.*[1 1 1])
    case 'regional_wbathy2'
        filename = [dirname, 'coast/pnw_coast_regional.mat'];
        load(filename);
        filename = [dirname, 'hybrid/pnw_combined_full.mat'];
        load(filename);
        contour(lon_topo,lat_topo,z_topo,[-4000:500:-500 -100],'color',0.6.*[1 1 1])
    otherwise
        error(['specify either detailed or regional or nothing'])
end

hold on;
plot(lon_coast, lat_coast, 'color', 'k', varargin{:})
clat = cosd(47.5);
set(gca,'dataaspectratio',[1 clat 1]);


